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Of 

The photon entropy production £ 7 is an extensive quantity and thus is fundamentally 
a volume integral of its local density <r 7 . Re-expression of £ 7 as a surface integral is effected 
through the divergence theorem and the surface photon entropy flux H. This step requires 
| the photons to propagate freely from a well-defined surface. Below this surface, the photons 
^ ' diffuse through matter instead. 
£ : 

In this case, of a star with a sharp surface of temperature T e g and no atmosphere, the 
paragraphs surrounding equations (3.10) and (A. 10) incorrectly present the surface integral 
form of £ 7 ("the total boundary entropy production, Ebound") as an additional contribution 
to £ 7 . In fact the boundary expression is simply the total S 7 within the enclosed volume, 
including all entropy production due to luminosity production and transport. The volume 
integrals (3.9) and (A. 10) are the photon entropy production up to, but not including, the 
surface; the surface integral or boundary term (3.10) is the photon entropy production 
including the surface, as well as the interior. 

Thus, the contribution of the surface alone to £ 7 is the difference of these two 
expressions: S surface = Sbound - S volume . The use of the volume integral (3.9) alone in the 
variational principle (3.11) for thermal structure continues to be justified by the boundary 
conditions that hold the radius R and T e g fixed. 
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The four equations of stellar structure are reformulated as two alternate 
pairs of variational principles. Different thermodynamic representations lead to 
the same hydromechanical equations, but the thermal equations require, not 
the entropy, but the temperature as the thermal field variable. Our treatment 
emphasizes the hydrostatic energy and the entropy production rate of luminosity 
produced and transported. The conceptual and calculational advantages of 
integral over differential formulations of stellar structure are discussed along 
with the difficulties in describing stellar chemical evolution by variational 
principles. 
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1. Differential Equations of Stellar Structure 

1.1. Static Stellar Structure as a Non-Equilibrium Steady State 

In thermostatic equilibrium (Reichl 1980; Callen 1985), general statistical properties 
radically simplify a macroscopic description. Non-equilibrium systems, however, generally 
require a microscopic kinetic theory. In special non-equilibrium systems, however, 
it is possible to leave the microscopic physics implicit in statistical averages, and to 
proceed with a nominally macroscopic dynamics. Stable stars evolve quasi-statically 
through well-defined spatial structures and temporal stages, so that a macroscopic non- 
equilibrium thermodynamics can be formulated. In this paper, we apply non-equilibrium 
thermodynamics to stellar structure, as distinct from stellar evolution. 

Macroscopic non-equilibrium thermodynamics began with the work of Rayleigh (1877) 
and Onsager (1931a, 1931b) (see also Casimir 1945); a complete theory was proposed by 
Prigogine and his followers which applies to an important but restricted class of systems 
(Prigogine 1945; Davies 1962; De Groot and Mazur 1962). Later developments have both 
extended (Donnelly et al. 1966; Glansdorff and Prigogine 1971) and rivaled this classic 
work (Tisza 1966; Truesdell 1969; Gyarmati 1970; Lavenda 1978). For stars, we may take 
a generalized Prigogine approach, assuming some type of local statistical equilibrium holds 
and intensive thermodynamics parameters are defined at least locally in space and time. 
This assumption validates a macroscopic approach. 

Non-equilibrium thermodynamics is typically described by conjugate pairs of variables: 
differences in intensive parameters, the thermodynamic forces, and macroscopic, extensive 
thermodynamic fluxes. Usually one relates forces to fluxes as cause to effect and assumes 
a linear or quasi-linear relation between the two. In this paper, we describe a quasi-static 
star of steady luminosity as an open system, receiving energy from nuclear sources which is 
ultimately radiated out of the stellar surface. 

For the mechanical and thermal structure of this non-equilibrium steady state (NESS), 
we present two different pairs of variational principles (equations (3.5) and either (3.11) or 
(4.2)). While the differential formulation of stellar structure integrates local quantities from 
point to point, either integral formulation directly starts with global properties, including, 
but not limited to, total mass, luminosity, and radius. Iterative application of a mesh 
approximation to the global variational integrals is analogous to the Henyey or relaxation 
differential method (Henyey, Vardya, and Bodenheimer 1965; Kippenhahn and Weigert 
1990). However, such a discretization approximation is not necessary, and continuous 
analytic approximations are also possible, in terms of variational parameters of global 
significance. 
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1.2. Formal Similarity Between Mechanical and Thermal Equations 

The four first-order differential equations of quasi-static stellar structure occur in two 
pairs: hydromechanical and thermal (Kippenhahn and Weigert 1990; Hansen and Kawaler 
1994). In the Euler representation and assuming spherical symmetry and conductive 
transport of luminosity, 



where the first horizontal pair is the mechanical (density-pressure) and the second the thermal 
(luminosity-temperature) equations. In radiative transport, the thermal conductivity k can 
be replaced by an equivalent radiative diffusion expression k — > 4acT 3 /3np, where k is the 
opacity of matter, ac/4 the Stefan-Boltzmann constant, and the Boltzmann constant ks is 
set to unity throughout. (Convective forms of luminosity transport are discussed below.) 
The variables are: r = distance from center; m(r) = cumulative mass from center to r, 
P(r) = P m {r) + .P 7 (r) = total pressure (matter and radiation), p(r) = matter density; 
l(r) = photon luminosity at radius r, T(r) = common matter-photon temperature, e(r) 
= luminosity production per unit mass, £j,(r)=neutrino luminosity per unit mass, G = 
Newton's gravitational constant. The luminosity is constant outside regions of luminosity 
production, where e = 0. 

Except in Section 5, we neglect the slow nuclear chemical evolution, so that the 
quasi-static stellar structure is in a NESS. Dropping the time derivatives, if the thermal 
conductivity were constant or oc T~ 2 , equations (1.1-2) would be self-adjoint and would 
admit a variational principle formulation. Because the conductivity or radiation transport 
coefficients generally do depend on dependent variables, the thermal equations are not 
self-adjoint. A self-consistent variational principle is nevertheless still possible (section 3). 

Except for this technical difference, the mechanical and thermal equations of quasi- 
static stellar structure are now formally symmetric: momentum transport and luminosity 
transport equations, luminosity conservation laws, transport coefficients and density are 
analogous. In this section, we exploit this structural similarity between mechanical and 
thermal equations, while still stressing the physical differences between the mechanical 
NESS and the steady photon entropy generation. 

In the static Euler representation, the seven dependent variables are: m, P, p, T, I, k, 
and e. There are four differential and three constitutive equations; seven equations in all. 
With four first-order differential equations, there are four boundary conditions which are 
usually taken as P(R) = (defining the surface r = R), m(0) = (no mass singularity at 
the center), m(R) = M (total mass of the star), and T(R) = T cS (surface temperature and 



—dP/dr = Gmp/r 2 
kdT/dr = l/Anr 2 , 




dm/dr = Anr 2 



P 



(1.1) 
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thus total luminosity). 

Because the included mass is conserved and accumulates monotonically with increasing 
radius, a point transformation to m as the independent variable is possible. In this Lagrange 
representation, the position r of each mass shell dm is now dependent, r = r{m). The 
structure equations now become: 

—dP '/ 'dm = Gmj '47rr 4 , dr / dm = 1 / 4nr 2 p , 
-kdT/dm = l/lQn 2 r i , dl/dm = e-e u . 

The boundary conditions in the Lagrange representation assume M as given. Thus, a useful 
set of boundary conditions is r(0) = 0, P{M) = 0, T(0) = T c , P(0) = P c . 



1.3. Differences Between Mechanical and Thermal Steady States 

The first and third equations in the Lagrange representation show one difference 
between matter and radiation and between mechanical and thermal structure. The star is 
defined by the presence of opaque matter, but, unlike m, luminosity / is still a dependent 
variable. Because it does not increase outside the core, / would be useless as an independent 
variable in the outside regions which transport, but do not produce, luminosity. 

A physically more important difference between matter and radiation arises because 
the matter is static so that matter thermodynamics enters only implicitly, through the 
equation of state P m = P m (p,T), where the matter internal energy acts as a potential for 
the pressure. The thermal structure is not static, but shows a stationary flow of released 
nuclear energy from the core to the surface. Luminosity transport is explicitly statistical 
and depends explicitly on thermodynamic averages: the local temperature T and the 
opacity K,(p,T). This non-equilibrium luminosity transport derives from non-equilibrium 
statistical mechanics (Appendix A). In radiative transport, this steady flow expresses the 
balance between the outward force of radiation and the resisting force of opacity (radiation 
friction): 

_dP 1 = KpJ_ , 13 , 
dr c 4irr 2 

Convection is more complex than radiative transport and usually described in mixing length 
theory by the mixing length A m i x as the transport parameter. The heat flux becomes: 

l/Anr 2 = (pc P T)X 2 mix ^ , (1.4) 

Ap 

for efficient convection, meaning that a convective cell loses little heat before it breaks up 
(Glansdorff and Prigogine 1971; Donnelly et al. 1966; Hansen and Kawaler 1994). (See 



- 5 - 



Appendix B for definitions.) Convection is efficient in two cases: The first is "slow but hot", 
where V is only slightly larger than V a d (quasi- adiabatic regime), and the heat transport 
is efficient because the heat capacity cp is large, but the cell velocity is slow. This regime 
occurs in the convective cores of stars with M > f.O8M . The second case is "warm but 
fast", where V is well above V a d, cp is not large, but the cell velocity is high. The outer 
convective layer of the Sun is in such a regime. 

Convection only comes into play as a luminosity transport mechanism if the 
Schwarzschild instability criterion holds. Convection, radiative transport, and conduction 
can operate simultaneously and in the same region of the star; in which case, the total 
luminosity is the sum of all three transport fluxes. In the case of efficient convective heat 
transport, virtually the entire luminosity flux is convective. 

The non-equilibrium steady-state obtains in quasi-static stages of stellar evolution, 
beginning with hydrogen and helium burning. Steady states, after the Main Sequence, burn 
carbon, neon, oxygen, or silicon, and are dominated by neutrino production and transport. 
The neutrinos escape essentially without matter interaction and without thermalizing. 
Neutrinos are produced by hot matter, but their energy loss is not thermostatic: their 
luminosity contributes to the material internal energy only by cooling the matter; they 
contribute nothing directly to stellar structure. The neutrino number and energy lost can 
only be inferred from the photon luminosity and the assumed matter properties and nuclear 
reactions (Bahcall f989; Kippenhahn and Weigert 1990; Hansen and Kawaler 1994; Arnett 
1996). 

2. Stellar Structure as a Non-Equilibrium Steady State 

No set of differential equations has a unique variational formulation (Douglas 1941). 
In this section and the next, we Legendre transform the global thermodynamic potentials, 
leaving the Euler-Lagrange equations invariant. The energy representation we start with 
enjoys a transparent physical interpretation in terms of hydrostatic energy and entropy 
production. In Section 4, we give entirely different variational principles in terms of local 
field variables, which nevertheless lead to the same mechanical and thermal equations (1.2). 

2.1. Thermodynamic Equilibria 

Thermodynamics distinguishes various types of equilibria or steady states (Gyarmati 
1970; Glansdorff and Prigogine 1971; Reichl 1980; Callen 1985). Global mechanical or 
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thermal equilibrium implies uniform pressure and temperature everywhere, or at least 
throughout a large, finite system. Stars are not in global equilibrium; their pressures and 
temperatures vary in space, but are still meaningful locally in small regions. If the length 
scales over which kinetic mechanisms maintain local mechanical or thermal equilibrium 
(LME or LTE) are much smaller than the scale heights of pressure or temperature: 



then a local pressure P(r) or temperature T(r) is well defined as a function of r. An 
equilibrium can be partial or complete among various constituents of the star. The LTE 
in a stellar interior is complete in that all matter species and radiation have a common 
temperature T m = T 7 = T. The pressure equilibrium is partial, in that the total pressure is 
the sum of the partial pressures of each component, P = P m + P 1 . The pressure decreases 
outwards because of the star's self-gravitation; the temperature decreases outwards wherever 
the total luminosity exceeds the neutrino luminosity. 

The local material internal energy density is determined by the matter pressure P m . 
In LTE, P m is an equilibrium state function of temperature and matter density, P m = 
P m (p,T), the matter equation of state. The photon equation of state is a function of T 
alone. In the luminosity equations, the transport coefficient (the opacity or conductivity) 
and the specific luminosity production, in LTE, are functions of local state variables 
T, p, and X h the local chemical composition (nuclear mass abundances Xj). The three 
constitutive functions depend on position only implicitly, through the local equilibrium 
variables T(r), p(r),Xi(r). Since we are generally not considering chemical evolution, we 
generally suppress the dependence on Xj. 

Both LME and LTE break down in the stellar atmosphere, where the matter becomes 
transparent and the radiation escapes almost unhindered. Here matter and radiation are 
not in thermal equilibrium. Pressure and gravity become unbalanced, as the star emits a 
stellar wind into space (Chandrasekhar 1950; Stix 1989) 

To distinguish these local thermodynamic equilibria from global pressure and 
temperature equilibria, the global state is referred to as a non- equilibrium steady state 
(NESS). Along the Main Sequence and any later, post-Main-Sequence steady states, 
NESS obtains, with steady radial luminosity flow. Various time scales are associated with 
these local equilibria and NESS's, which must be achieved and maintained by different 
mechanisms. In the Sun, LME and LTE are achieved in about 10~ 12 sec by local kinetic 
mechanisms. The global mechanical and thermal NESS's are established in about 30 
min and 10 7 yrs, respectively. These last two time scales characterize helioseismological 
disturbances and non-static, macroscopic heat flows (Kelvin-Helmholtz time). Except in 
dynamic stages, these time scales are much shorter than the local and global chemical 





(2.1) 
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evolution times 10 7 8 and 10 10 yrs associated with nuclear transmutations (Bahcall 1989; 
Kippenhahn and Weigert 1990; Hansen and Kawaler 1994; Arnett 1996). 

2.2. Different Thermodynamic Representations 

Like any thermodynamic system, a star can be described by different global 
thermodynamic potentials. The original representation is the energy representation E = 
E(V, S, Ni), as a function of volume V, entropy S, and species numbers iVj. Because the 
thermodynamic state varies spatially, the extensive state variables must be locally recast, 
as either spatial (per unit volume) or specific (per unit mass) densities (Chandrasekhar 
1939; Kippenhahn and Weigert 1990). (We follow astrophysical custom by using the latter, 
unless otherwise noted.) Extensive quantities are then mass integrals over the specific 
densities. The specific energy e is the sum of the specific internal energy u and the specific 
gravitational energy —Gm/r. The specific volume v p is l/p. In all but the hottest stars, the 
total entropy S — S m + S 1 is dominated by matter, with specific entropy s m . Nevertheless, 
the small radiation entropy flux is responsible for the stellar luminosity and cannot be 
ignored. 

In the stellar interior, where a NESS obtains, the specific entropy s = s m + s 7 is 
stationary, so that E = E(V, S) and e(p, s m + s 7 ) are, respectively, global and specific 
thermodynamic potentials. The radiation entropy density alone is determined by the local 
ratio nl/m. If the thermal structure of the star is prescribed (e.g., isentropic or isothermal), 
a barotropic relation P = P(p) holds so that the hydromechanical equations are closed. 
Nuclear burning in the stellar core leads to secular chemical changes (which we ignore 
in stellar structure), and to the rapid thermalization of fusion products, which steadily 
produces radiation entropy flux out of the star. In the steady state, the star's entropy is 
constant and cannot be used as a global thermodynamic state function. In order to treat 
matter and radiation as symmetrically as possible, we use the temperature T, which is 
common for the matter and radiation in the stellar interior. The appropriate generalization 
for the stellar atmosphere is straightforward. (See Appendix A, subsection 3.) 

The thermodynamic potential for the new variables (V, T) is the Helmholtz free energy 

F(V, T) — E(V, S) — TS (2.2) 

or specific Helmholtz free energy f(p,T). (We might also go over to completely intensive 
state variables (P, T) and the Gibbs free energy G(P, T) = F + PV as the thermodynamic 
potential. The choice of thermodynamic potentials is a matter of convenience; we prefer the 
Helmholtz free energy.) 
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3. Global Thermostatic Potentials and Entropy Production 

In this section, we denote functional variations by 5, local thermodynamic state and 
spatial changes by ordinary differentials, so that the spatial differential dF = VF • dr. 



3.1. Mechanical Steady State 

In considering only the mechanical NESS, we can use any thermodynamic 
representation. We can use the total hydrostatic energy (Lamb f 945; Chiu 1968; Rosenbluth 
et al. 1973; Hansen and Kawaler 1994): 

r-M 

E(V,S)= dme(p(r),s(r)) , e = u-Gm/r , (3.1) 
J o 

or the Helmholtz free energy 

r-M 

F(V,T) = dm f(p(r),T(r)) , f = e-Ts , (3.2) 
J o 

as the thermodynamic potentials. The specific internal energy and free energy have 
variations 

5e = -P5(l/p)+Gm5r/r 2 + T5s , (3.3) 
5f = -P5(l/p) + Gm5r/r 2 - s5T , (3.4) 
where 5(1/ p) = d(4nr 2 5r)/dm. 

We could also use the global enthalpy 

H(P, S) = E(V, S)+ J dV P , dH = VdP + TdS , 

or the Gibbs free energy 

G(P, T) = E(V, S) + J dV P - J dS T , dG = VdP + SdT , 



as global potentials. The hydrostatic equation (1.2) is obtained by freezing the thermal 
structure, taking adiabatic variations in E or H or isothermal variations in F or G. For 
either mechanical variation, 5P/ p = 5h or 5g. In all cases, the functional variation does not 
assume that the frozen thermal structure is in the NESS. But the short mechanical time 
scale compared to the thermal time scale assures that LTE, although not necessary, is a 
sufficient condition for LME. "Isothermal" here refers not to spatially constant temperature, 
but unvaried temperature profile; similarly, "adiabatic" refers to unvaried specific entropy 
profile. 
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The simplest example of the resulting Euler-Lagrange equations arises from applying 
the variation (3.3) to (3.1), setting 5s = (Chiu 1968). Using 5(1/ p) = d(Arrr 2 5r)/dm and 
integrating the first term in (3.3) by parts, one obtains 



ignoring the boundary condition terms. Setting the integrand here to zero for arbitrary 5r 
reproduces the hydrostatic NESS, the first equation of (1.2). Alternatively, applying the 
variation (3.4) to (3.2) isolates the mechanical structure in the isothermal limit 5T = 0. 



In our NESS, the entropy density of matter and of radiation are each stationary and 
very disparate in magnitude. Because matter and radiation share a common temperature, 
we use T as the intensive thermodynamic variable. Instead of the entropy, our variational 
principle for thermal NESS minimizes the entropy production rate E > (Prigogine 1945; 
Davies 1962; De Groot and Mazur 1962; Gyarmati 1970; Donnelly et al. 1966; Glansdorff 
and Prigogine 1971). A simple minimum entropy production variational principle applies 
to systems close to global equilibrium in the linear regime and to a special class of systems, 
those in the quasi-linear regime. In the linear regime, transport coefficients are constant. 
In the quasi-linear regime, LTE still obtains, but the constitutive functions depend only on 
the local state variables and not on their derivatives. Because the stellar NESS is in the 
nonlinear regime, the thermal variational principle depends on two temperatures, T, T*, 
the first of which is dynamical, the second referring to the background matter. 

The entropy production rate is bilinear, a product of thermodynamic forces and 
fluxes. For temperature variations, the thermodynamic force arises, not from T, but from 
differences of 1/T, or from the gradient of l/T. In the linear or quasi-linear regime, the 
luminosity flux is linear in V(l/T). This thermal gradient drives the flow of luminosity from 
point to point, for any of the three standard heat transport mechanisms. 

The temperature also occurs in other parts of the entropy production rate, but only 
as a local state variable having nothing to do with heat transport. We must distinguish 
this temperature, T*, from the temperature T which is subject to functional variations 
(Donnelly et al. 1966; Glansdorff and Prigogine 1971). With this distinction, we can write 
the entropy production of conductive transport: 




(3.5) 



3.2. Entropy Production in the Thermal Steady State 




(3.6) 
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and of radiative diffusion as: 

r 1 AnrT 5 

s 7diff = y ^ 2 ( ^7 } * [v(i/t)]2 ' (3 - 7) 

with details in Appendix A (Chandrasekhar f950; Essex f984a, 1984b). The entropy 
production of convective transport is presented in Appendix B (B.2). Any transport entropy 
production is a measure of the efficiency of these different mechanisms for luminosity 
transport, the price paid for transporting the luminosity produced in the core of the star to 
the surface. 

The entropy production due to nuclear burning in the core is: 

S fusion = [epl/T . (3.8) 

The bulk radiation entropy production rate is the sum of the luminosity transport and 
production terms (3.6-8): 

^7 ^7 difT ^fusion /q q\ 

= JdV (i/2)[4acT 5 /(3 K p)]4V(l/T)] 2 + M,/T , l ' J 
in case of radiative transport. 

Where the radiation outstreams from the photosphere, at temperature T e g, 

Sbound = ^(00/4)7^(47^) (3.10) 

is obtained by integrating the radiation entropy flux over the surface of the star. Because 
the mechanical NESS is reached so much faster than the thermal NESS, the mechanical 
structure must be the hydrostatic NESS corresponding to a given thermal distribution, 
stellar radius R and surface temperature T e g (Sieniutycz and Berry 1989, 1991, 1992, 1993; 
Sieniutycz and Salamon 1990; Sieniutycz 1994). This boundary entropy production is held 
fixed by the boundary conditions and plays no dynamical role in the thermal variational 
principle. 



3.3. Minimum Entropy Production 

Our thermal variational principle requires minimizing the entropy production by 
varying T — > T + ST, but not T*, holding M and T* fixed (Donnelly et al. 1966). Only after 
the Euler-Lagrange equation is obtained, is T* set equal to T. For each of the transport 
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entropy productions (3.6, 3.7, B.2), the standard heat transport equations are obtained in 
this way. For example, the radiative diffusion case is derived from (3.9): 



c r 4-nrT 5 n 
,)S 7 = jdV [[^— UV(1/T)-W(1/T)] - [epU5T/T 2 )} = . (3.11) 



Integrating the first term by parts and dropping the boundary terms, we set T* to T and 
the entire integrand to zero. For arbitrary variations ST, the result is radiative diffusion: 

d A6nr 2 acT 3 dT 



, s o ) • ^r 2 ep . (3.12) 

dr v 6np dr ' 

For simplicity, only the radial dependence has been kept, and the neutrino luminosity e v 
has been ignored. This same result is obtained by combining the thermal-luminosity pair of 
equations in (1.1). 

In practice, an ansatz is made for Tjj) and for T(r), with only the latter containing 
variational parameters to be determined by minimizing the variational integral. The 
procedure is iterative: after each step T*(r) is set equal to the T(r) obtained at the previous 
step. One may use a global mesh, analogous to that in the differential Henyey method, or 
an analytic form with a finite or infinite number of adjustable and physically suggestive 
parameters. The procedure is identical to the background field (BF) or self-consistent field 
(SCF) method used in quantum many-body and field systems and can lead to an analytic 
approximation for the stellar structure in terms of global parameters. These may be related 
to global properties and boundary conditions of the star. 

Because no entropy production is associated with the thermalization or diffusion of 
neutrinos, the only neutrino entropy production is by averaging over neutrino energies 
produced by the thermalized reactant matter. The exact expression for neutrino entropy 
production varies by reaction, but for high temperatures, its dimensional order of magnitude 
is very approximately: 

H v ~ N v , 

where N v is the total neutrino production rate. (A better estimate is given in Appendix C, 
C.10) The neutrino contribution to entropy production is usually much smaller than the 
photonic. But in advanced stages of stellar evolution, the temperature is high enough to 
bring weak interactions significantly into play, letting neutrino entropy production rival 
that of heat and radiative diffusion mechanisms. 

In supernova explosions, neutrinos do play the same role as photons in ordinary stars: 
they interact with the ambient stellar matter, are thermalized, and are then emitted from 
the surface of the neutrinosphere at a definite blackbody temperature (Bahcall 1989; Arnett 
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1996). The temperature in the innermost region reaches a very high 30-50 MeV, so that 
neutrino luminosity transport is efficient and neutrino entropy production remains fairly 
small, below that of other energy transport mechanisms. 

4. Variational Principles with Local Specific Potentials 

In the preceding section, we presented variational principles for the hydrostatic NESS 
(3.1-5) and for steady-state luminosity transport (3.6-8). For mechanical NESS alone, 
any global thermodynamic potential can be extremized, but for the thermal steady state, 
only the global potentials, G(P,T) or F(V,T), and a close relative, S(T), could be used. 
Independently of thermodynamics, we know that the same differential equations admit of 
many different variational integral forms. Because the specific potentials, potentials per 
unit mass, vary in space and also represent the local thermodynamic state just as well as 
the variables (V, P, T, S), we may take the specific potentials as field variables and find 
other variational integrals to serve as global actions. 

In the mechanical equations, we can switch from the density p and pressure P as the 
intensive state variables to either the specific enthalpy 

h(P, s) = e(p, s) + Pp , dh = p dP + Tds (4.1a) 

or the specific Gibbs free energy 

g{P,T) = e{p,s)+Pp-Ts , dg = p dP + sdT , (4.1b) 

as the new local field variables. The mechanical NESS is obtained by either varying 

^ad = JdV h(^)(V/>) 2 + P(h)] (4.2a) 

adiabatically or 

2L = J dV ^^G^gf + P(g)] , (4.2b) 

isothermally, so that bPj p = either Sh or Sg. 

For the thermal part, we define a heat or radiation potential density 8, where d9 = 
Jq ■ dr and dr is the distance increment through which the luminosity is transported. In 
conductive or radiative transport, d6 = —kdT or — (c/ Kp)dP 1 . In convective transport, 



d6 = (pc P T)(A mix /Ap)[V - V ad ](u> dr) 



(4.3) 
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where w = upward speed of the convective cell, assuming V > V a d (see Appendix C). The 
luminosity source is included by a source potential 11(0), where dU = pedO. The thermal 
variational integral is then: 

J = JdV[-(±-)(Wf-n(6)] , (4.4) 

where r is the independent variable and 9 the field. Variation of (4.4), keeping the 
mechanical structure fixed, yields the thermal transport equations for all three kinds of 
transport. These actions X, J do not have the simple thermodynamic interpretation that 
the minimal entropy production principle enjoys. 



5. Variational Thermohydrodynamics and Chemical Evolution 

A full statement of stellar evolution using entropy production and variational principles 
is beyond the scope of this paper, and only a sketch is presented in this section. Evolution 
refers only to nuclear chemical evolution, not to faster hydrodynamic or thermal changes 
with fixed nuclear abundances (Kippenhahn and Weigert 1990; Arnett 1996). The latter 
two processes are only time-dependent versions of the hydrostatic and thermal NESS. Such 
dynamical but non-evolutionary behavior can be treated with time-dependent extensions of 
the static variational principles stated in section 3. The hydrodynamic part is derived by 
minimizing the action (Lamb 1945; Sieniutycz and Salamon 1990; Sieniutycz 1994): 

/^/rfm[^-i(V$) 2 + e(p, S )] , (5.1) 

where <£> is the velocity potential, v = — V$ = the macroscopic velocity of fluid flow. This 
action principle is adiabatic and holds only for irrotational fluid flow. The isothermal 
analogue of (5.1) can be obtained by substituting f(p,T) for e(p,s). (A hydromechanical 
small-oscillation variational principle has been derived for asteroseismology by Backus 
and Gilbert 1967). The non-steady thermal behavior can be derived by minimizing the 
variational integral (Donnelly et al. 1996; Glansdorff and Prigogine 1971): 

J dt J dV (X C ond + 0" 7 diff + or COIlv + ^ + [pc v T%^(^j , (5.2) 

where cy is the specific heat capacity at constant volume, assuming the mechanical part 
to be instantaneously in the NESS (Sieniutycz and Salamon 1990; Sieniutycz 1994). 
Equation (5.2) is not a minimum entropy production principle, because it represents a 
NESS only if dT/dt = 0. 
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The two variational principles (5.1,2) can be be used to treat stellar oscillations, with 
or without linearization. In general, these principles encompass only the propagation of 
seismic waves, not their driving forces, unless the latter are explicitly included in the 
integrals. 



5.1. Gravitational Settling in Chemical Steady State 

Stellar evolution involves two fundamental changes: gravitational settling of heavier 
nuclei arising from spatially inhomogeneous nuclear fusion reactions and thermonuclear 
transmutation of elements. In the Sun, for example, the associated global time scales are 
6 x 10 13 and 10 10 yrs, respectively. Both evolution processes can be described in terms 
of entropy production, but only element diffusion is close to local chemical equilibrium 
(Kippenhahn and Weigert 1990; Bahcall and Pinsonneault 1992). As the mechanical and 
thermal timescales are usually much shorter than the chemical timescale, the star is at each 
instant in evolution in both hydrostatic and thermal NESS. But the boundary conditions 
for these NESS's change with chemical evolution, as the star's mass, radius, and luminosity 
change. 

Chemical diffusion by gravitational settling can be dynamically formulated in terms of 
a minimum entropy production S nuc diff (Donnelly et al. 1966; Glansdorff and Prigogine 
1971): 



J nuc diff 



JdV 



+ 



E y |[^(p,T)],V(™)v(^) 

E^t(p,t)],v(™)v(i) 



(5.3) 



where the pi(r,t) are chemical potentials of the nuclear species i, and the T>ij(p,T) are the 
chemical diffusion coefficients, including cross terms between different species. The mixed 
thermo diffusion effect is included with the V iT terms. By Onsager's theorem, = T>^ 
(Onsager 1931a, 1931b). Each chemical potential is defined relative to some reference 
potential, here taken to be that of hydrogen (H). This form of chemical diffusion assumes a 
quasi-linear relation between the element flux Jj and the chemical gradients V//j : 

^ = -J2^MT)V(^^)-V tT (p,T)V(^) , (5.4) 

j 

in analogy with heat conduction. The diffusion represented by (5.3) is therefore in a NESS, 
as equation (5.4) is time-independent. 

Chemical diffusion could also be formulated in terms of the extensive variables N iy the 
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species densities rij, or the nuclear mass abundances JQ, by use of the grand potential: 

A(V,T,fMi) = F-JdVJ2iVini , , g g v 

Hi = -p(da/dfii) , 

where a(p,T,pi) is the specific grand potential. 



5.2. Thermonuclear Burning as Nonlinear Reactions 

The entropy production of nuclear fusion is composed of a radiation part E 7 , already 
discussed in section 3, and a matter part E matte r, hitherto ignored, because it is significant 
only in late stages of stellar evolution. Starting with the Gibbs form of the entropy 
increment (De Groot and Mazur 1962; Gyarmati 1970; Glansdorff and Prigogine 1971), the 
canonical expression for £ of matter undergoing chemical/nuclear reactions is: 

^matter = J dV -JT ■ ( 5 - 6 ) 

where n,i are the number densities of each species i. The photon chemical potential p 1 is 
zero, because photon number is not conserved. The neutrino chemical potential p v is zero 
as long as neutrinos have no significant density. 

A more transparent form can be obtained by defining common reaction rate densities 
£,ot = ni/v a ,i, where the subscript a labels the reaction. The v a i are the stoichiometric 
coefficients of reaction a : 

i/ a ,iNi + v a ,2N 2 + ... -> -is a ,j+iNj - v a j +2 N j+2 + ... , 

for all participating reactants i = and products i — j + 1.... The affinity A a of reaction 
a : 

Aa = Yl "acifr ( 5 - 7 ) 
i 

is a measure of how far a given reaction a is from chemical equilibrium. For a dead star 
(one whose nuclear reactions have gone to completion), A a = for all a. The alternate 
form for E matter is then: 

S mato = / dV J2 ^ , (5-8) 

a canonical bilinear form in the forces A a /T and fluxes £ Q . 

Because A a ~ MeV is higher than the ambient temperature T, the thermodynamic 
forces driving nuclear burning are large. This makes stellar thermonuclear reactions highly 
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nonlinear, even locally (Reichl 1980). Because no general linear or quasi-linear relation 
holds between the fluxes and the forces in the matter contribution to E, without additional 
knowledge or details, we can go no further in expressing the entropy production of matter 
due to nuclear fusion than the kinematical expressions (5.6, 8). Without a quadratic (or 
any other expression) for E mat t er in terms of the forces A a /T there is no automatic principle 
of minimum entropy production (Prigogine 1945; Lavenda 1978). 

A variational principle of the quadratic type requires a quasi-linear relation between 
forces and fluxes. A non-quadratic minimum entropy production principle for matter 
fusion may still be possible, if the chemical evolution is in a NESS, with the fluxes £ Q 
constant in time. This NESS approximation holds as long as the star remains on one 
branch of hydrostatic thermonuclear burning (H, He, C, Ne burning), or on one of the late, 
hydrodynamic stages of chemical evolution (O, Si burning) (Arnett 1996). Furthermore, 
each burning stage constitutes a different NESS. It is possible that minimum E mattcr holds 
on each burning stage, but that the bifurcation from one burning stage to the next is a 
discontinuous phase transition similar to the switch from non-convective to convective heat 
transport (Appendix B). So long as some fuel remains locally from one stage, that burning 
continues; the minimum entropy production conjecture would imply that only when its 
fuel runs out is the next burning stage preferred on entropy production grounds. Similar 
considerations may also hold for the evolution of protostars. 

If the matter fusion is not representable as a NESS, another, time-dependent variational 
principle, analogous to (5.2), may still be possible, lacking a simple interpretation as 
minimum entropy production. 

6. Summary and Outlook 

Variational principles are most useful in expressing general properties of static or 
dynamical systems, such as exact or approximate symmetries, and in suggesting generalized 
coordinates. Integral principles also open alternative paths to solution via an iteration 
of approximate solutions. Under some restrictions, such methods converge and provide a 
useful replacement for numerical differential methods (Donnelly et al. 1966). 

We have presented two different pairs of variational principles to replace the standard 
four equations of a hydrostatic NESS and steady luminosity flow. Two of these, the thermal 
variational principles (3.6-8) and (4.4) for NESS luminosity transport, are the principal 
original results of this paper. These represent a global alternative to numerical integration 
of the four differential equations of stellar structure. Practically, these thermal variational 
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principles suggest calculational procedures analogous to either the relaxation technique or 
the Rayleigh-Ritz method (Donnelly et al. 1966; Glansdorff and Prigogine 1971). Examples 
of global analytic approximations will be published elsewhere. 

As seen in section 5, the static variational integrals of sections 3 and 4 can be extended 
to the time-dependent hydrodynamic and thermal non-steady states. 

We have been able to recast stellar evolution only partially into integral form, and 
a dynamical principle is lacking. A complete integral reformulation of stellar theory 
would encompass these outstanding issues: the transmutation of elements, the emission of 
neutrinos, and bifurcations and instabilities during the multiple burning stages. 

We are indebted to Christopher Essex (Univ. W. Ontario) for many helpful discussions 
and suggestions. This research was supported at the Institute for Theoretical Physics, 
University of California, Santa Barbara, by the National Science Foundation under Grant 
No. PHY89-04035; by the University of Florida, Institute for Fundamental Theory of 
the Department of Physics; by the Aspen Center for Physics; by the Telluride Summer 
Research Center and Telluride Academy; and by the Department of Energy under Grant 
Nos. DE-FG05-86-ER40272 (Florida) and DE-AC02-76-ERO-3071 (Penn). We thank the 
ITP, the Aspen Center, and the TSRC for their hospitality. 

Appendix A: Conductive and Radiative Entropy Production 

Let us define notation and simple concepts first. LTE is assumed to be valid in the 
strictly local limit, with small deviations over small distances. Thus matter and radiation 
are in the maximal entropy state locally, with non-zero second-order thermodynamic 
fluctuations about the LTE. Purely local state variables and functions carry an asterisk (*) 
subscript. These are not subject to the non-equilibrium functional variations (Donnelly 
et al. 1966; Glansdorff and Prigogine 1971; Chandrasekhar 1939 and 1950; Essex 1984a, 
1984b; Holden and Essex 1996). 

1. Heat Conduction 

The familiar case of heat conduction best introduces the entropy production E. This 
function's spatial density a takes on, for matter alone, a simple bilinear form derived 
from the Gibbs equilibrium formula for the entropy increment dS = dQ/T + mechanical, 
chemical,... etc. terms (Callen 1985; Balian 1992). The entropy produced by a heat current 



Jq flowing through a temperature field T in a volume V with a surface area A is: 

S cond = JdA-(3 Q /T) 
= jWV-(J Q /T) 

In a steady state, V • Jq = q, where q is the heat source density within the volume. On the 
other hand, Jq in the quasi-linear case is linear in an externally given inverse-temperature 
gradient, with: 

J Q = (A;T 2 )V(1/T) , (A.2) 

where k(p, T) is the local thermal conductivity, an integral over the microscopic momentum 
space: 

kT 2 = (47r/3) / dpp 2 f (r,p)[v 2 (pe) 2 /w(p)} , (A.3) 
Jo 

where v and p are the atomic speed and momentum, pe the gas internal energy density, /o 
the reduced 1-particle phase space distribution at LTE, and w(p) the collision time (Balian 
1992). Thus: 

<Tco„d = (kT 2 )*[V(l/T)] 2 + q*/T (A.4) 
is the entropy density production rate by the transport and production of heat. 

In (A.4) the current Jq must be evaluated assuming a given external T gradient. In 
reality, this gradient is produced by the current itself. Thermo dynamically, the form (A.4) 
assumes an externally-imposed first-order deviation in entropy from LTE, while in fact, 
the deviation for the transport term is a second-order fluctuation (Donnelly et al. 1966; 
Glansdorff and Prigogine 1971). (LTE makes the first-order fluctuations of entropy zero.) 
Correction of this problem leads to an extra factor of 1/2 in the first (transport) term. Note 
that this term is quadratic in the gradient. The second (production) term receives no extra 
factor of (1/2), as it really is a first-order deviation in entropy: the heat production q is an 
externally given function. Fourier's equation for heat conduction is obtained by varying T 
in 

Scond = JdV (l/2)(kT 2 )JV(l/T)] + q./T\ , (A.5) 

holding (kT 2 )* and fixed. Afterwards, T* is self-consistently set equal to T. The heat flux 
is 

J Q = (//47rr 2 )f = -WT . 

2. Radiative Diffusion 

Radiative transport by photon diffusion is formally similar to heat conduction by 
matter-matter collisions (Chandrasekhar 1950; Essex 1984a, 1984b; see also: Planck 1913; 



(A.I) 
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Rosen 1954; Kroll 1967; Lallemand and Martinet 1979). The role of k is taken by an 
expression involving the opacity k of the matter to photon travel. We should then expect 
k to involve an integral over the photon phase space. The evaluation of S 7 for LTE with a 
small gradient begins with the generalized bilinear form for photons at angular frequency uo 
passing through and interacting with matter at temperature T : 



[OO f + l 

2tt jf dcu J i di J u [l/T u - l/T] , (A.6) 



>j difr 



where £ = photon local direction cosine ( not to be confused with the reaction rate densities 
£ a of section 5), is the differential radiation luminosity density out of equilibrium: 

Jw = K w p\B w - I u ] , (A. 7) 

with K/jj the frequency-specific opacity of matter, B u the Planck function (blackbody 
differential radiation energy flux), and I w the true energy flux of photons. In the spherical 
diffusion approximation, I u = B u - (£ / n UJ p)(dB u} / dr) with the gradient term small except 
very near the stellar surface. T u is the effective brightness temperature for any l w and varies 
with uj : 

l/T w = (l/Mln[^ ] - + l] • (A.8) 
In the Rosseland mean opacity: 

l/« = (^°° duo [1/k w ] dB u /dT)l(j™ duo dB w /dT) , (A.9) 

the denominator has the value acT 3 . 

The entropy production is a generalization of the bilinear Gibbs formula, with 
changes in l/T playing the role of the thermodynamic force (which is functionally varied) 
multiplying some flux, as outlined in Section 3. 



3. Source and Boundary Terms 

Nuclear burning produces entropy by the production of high-energy photons and the 
thermalization of fusion products. The radiation/kinetic energy is produced in a fusion 
reaction by thermalized matter, a tiny, positive contribution to radiation entropy, as the 
original matter reactants are thermalized. This original photon/kinetic energy is absorbed 
upon thermalization, a negative contribution to entropy. Both the matter kinetic energy 
and radiation are then thermalized to the ambient temperature of the core, a large and 
positive contribution. The first two contributions are negligible compared to the third, 
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being suppressed by the ratio T/T , where T is the brightness temperature of the original 
photons, in the range 0.1-5 MeV, well above the typical stellar core temperature. These 
contribution are small but non-negligible for older stars with higher core temperatures. 

For any kind of radiation transport, there is a constant entropy production from the 
release of radiation into empty space. The local radiation entropy flux H has magnitude 
(4/3)(ac/4)T c 3 ff on a stellar surface with temperature T e g. Multiplying this expression by the 
surface area 4tiR 2 gives the total boundary entropy production, Sb OU nd- As in conduction, 
the transport term arises from second-order fluctuations in the LTE entropy; thus the factor 
of 1/2 again. The luminosity source term ep is again a true external first-order deviation 
from LTE, so the complete entropy production is: 



in the radiative diffusion case, where T* and M are not varied. 

A more realistic procedure would be to construct the E 7 for the stellar atmosphere. 
Multiple matter temperatures and a variable radiation brightness temperature must then 
be introduced. We do not discuss stellar atmospheres in further detail. 

A sense of the relative sizes of these entropy production rates may be gotten from 
the example of the Sun: Conductive heat transport is small; the core and radiative zone 
transport luminosity by radiative diffusion. We estimate the various entropy production 
rates, using the boundary contribution E bound ~ 8xl0 29 erg/°K sec as a benchmark. The 
radiation entropy production in the core is much smaller, because of the much higher 
temperature T cS ~ 5500°K versus T c ~ 15 x 10 6o K : E fusion ~ (0.0007)E bound . The entropy 
production rate due to transport of the luminosity from the core to the surface is E di g — 
(0.0002)E bound . This contribution is also much smaller than the boundary term, but not 
much smaller than the core term. Essentially all of the radiation entropy production in the 
Sun arises from the release of radiation into empty space at the surface. But, since this 
term arises from boundary conditions, it is not varied in deriving the transport equation. 



Unlike conductive and radiative transport, convection is a macroscopic process, 
involving bulk motion of matter. Furthermore, real convection is complicated, with many 
length scales and the possibility of significant turbulence effects. No satisfactory theory 
of convection exists, but there are useful models that capture the essentials of the heat 
transport, enough for consideration in the NESS structure of MS stars (Chandrasekhar 



E 



'7 




(A.10) 
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1939; Lamb 1945; Donnelly et al. 1966; Glansdorff and Prigogine 1971; Stix 1989; Sieniutycz 
and Berry 1989, 1991, 1992, 1993; Hansen and Kawaler 1994). 

1. Schwarzschild Criterion and Mixing Length Theory 

Schwarzschild's well-known picture begins by idealizing the convective bulk motion of 
matter cells slightly hotter than their surroundings and heated from below in a gravitational 
field by a temperature gradient. In the basic mixing length theory (MLT), such a cell, 
with linear size A mix , moves upwards a distance A mix before breaking up and merging 
with the surrounding gas or fluid, releasing its heat in the process. (In real convection, 
there are many mixing scales, not one.) The cell floats upwards by bouyancy, its internal 
temperature T' being higher and density p' lower than its surroundings. In MLT (Boussinesq 
approximation), the effects of sound and shock waves are ignored; the pressures internal 
and external to the cell are assumed equal; and the T and p variations between the inside 
and outside of the cell are small: (p — p')/p and (T" — T)/T <C 1. In stellar applications, 
in addition, A m i x is assumed much smaller than the stellar radius R (Hansen and Kawaler 
1994). 

In order for convection to occur at all, the temperature contrast between the bottom 
of a mixing length, where the cell is heated, and the top must be large enough to overcome 
gravity by creating a bouyancy force. In the continuum, this leads to Schwarzschild's 
criterion for convection to occur: 

V>V ad , (B.l) 

where V = — d\nT/d\nP, the temperature T(r) and pressure P(r) are the actual profiles 
of a given star. V a d is the same as V, evaluated for the same star, but as if the equation of 
state were adiabatic, which assumes that the cell exchanges heat at most very slowly with 
its surroundings before it breaks up. In the convective regime, the temperature gradient in 
the star outwards is steep enough that relatively hotter cells are bouyant. The simplicity 
of this condition hides an important complication: the presence or absence of convection 
itself changes the T(r) and P(r) profiles. In practice, the problem is solved iteratively and 
self-consistently. For quasi- adiabatic convection V is only slightly larger than V a d, and 
the details of this "slow but hot" convection are unimportant. For V well above V a d, the 
specific assumptions of the convection model become important, including turbulence in 
stellar conditions with very low viscosity (very high Reynolds number) and smallness of 
Amix compared to the pressure scale height Ap. Such convective zones are then "warm but 
fast"; but, typically being in the outer parts of stars, they are almost irrelevant to the 
stellar cores. 
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2. Entropy Production: Heat Loss and Bouyancy 

As a convective cell rises, it can lose part of its excess heat by thermal 
conduction/radiative diffusion (the two can be treated on an equal footing by using 
the equivalence of Appendix A) and by viscosity The total entropy production of 
convection E conv is a sum of three terms: the first due to heat loss (positive), the second 
to viscosity (positive), and the third to bouyancy (negative). The last is negative because 
bouyancy converts heat to the work of raising the cell in the local gravitational field "g" . 



""conv "heat loss "vise "b 



ouy 



2 



\2 



"heat loss = [(pC P ) 2 /(2/c)]|wA mix /ApJ (V - V ad 
"vise = r?U>7( T ^mix) , 

"bou y = -[2(pg 7 p)(wA mix ) 2 /(r ?T TAp)] (V- V ad ) , 
where w is the upward velocity of the cell. In MLT, w is solved for self-consistently: 



(B.2) 



w = \ mix y (ipg/ Ap) (V - V ad ) . 

The kinematic viscosity is rj, the thermal diffusivity r\ T = k/ pep, and cp is the specific heat 
capacity of matter at constant pressure. The isobaric exponent 7p = — (din p/d In T)p. The 
temperature contrast of the cell (T") with its surroundings (T) is 

The convection is efficient if the cell's cooling time X^ix/Vr is long compared to its travel 
time Xp/w. In stars, the viscosity is still smaller and negligible in MLT, i.e., r] <C tjt- 

The inverse temperature difference that acts as the thermodynamic force for heat 
transfer is just the expression (B.3). The relative temperature contrast, T times 
equation (B.3), must be small for quasi-linear thermodynamics to be valid. The range of 
its validity is identical to that of MLT. In the quasi-adiabatic or "slow but hot" regime, 
V — Vad is positive but so tiny that the relative temperature contrast is small. In the "warm 
but fast" regime, V — V ad is positive and not tiny, and the relative temperature contrast is 
not small. Quasi-linearity and MLT are then not valid, although MLT is commonly used 
anyway, for the lack of a better but still simple model. The specific details of convection 
also matter in this regime, and the simplicity of the quasi-linear situation vanishes. 

Convection and convective instability are illuminated by comparing the competing 
entropy production rates for the same cell in two regimes, convective (heat loss + bouyancy) 
and non-convective (heat loss alone) (Figure 1) (Donnelly et al. 1966; Glansdorff and 



-23- 



Prigogine 1971). Let us start by taking a convective motion with an arbitrary cell velocity 
w. Without the bouyancy part, the entropy production due to heat conducted out of the 
cell falls and then rises quadratically, reaching zero at V = V a d (Figure 1: solid line). (The 
cell velocity is actually zero in this case and must be replaced by an equivalent expression 
for pure conduction.) If we put the bouyancy part back in, the entropy production curve 
is modified by a negative term linear in V — V a d and becomes distorted (Figure 1: dashed 
line). The two curves cross at zero, for V = V a d- The branch with lower entropy production 
for V < V a d is the heat loss alone. But for V > V a d, the heat loss + bouyancy branch 
has the lower entropy production rate and is favored. In this regime, the cell becomes 
bouyant, and convection begins. The analogy with first-order equilibrium phase transitions 
is evident, with the entropy production rate taking the place of the free energy. 

In the case of convection, there is no variational calculation for the thermal structure, 
only a comparison of two discrete branches of S. The cell velocity w must be computed 
from hydrodynamic considerations, a topic beyond the scope of this paper. 



Because untrapped neutrinos are not in LTE, they do not contribute directly to stellar 
structure. While the total lepton number (charged leptons plus neutrinos) is conserved, 
neutrino number is generally not conserved. After their production in the hot stellar 
core, in quasi-stellar stars, neutrinos suffer negligible weak interactions and change in 
neutrino number. They therefore transport both neutrino luminosity and number, while 
electromagnetic radiation transports only photon luminosity. Neutrinos contribute to 
quasi-static stellar structure only indirectly by cooling of matter. 



In analogy with the specific photon/heat luminosity production rate, a specific neutrino 
rate, e u , is defined. As for photons, this quantity is a sum over all neutrino-producing 
reactions, 



where is the neutrino differential production rate density and we label the neutrino 
phase space by neutrino energy E, not frequency. The total specific luminosity is the 
function e. But unlike photons, neutrinos are not thermalized after they are created, leaving 
the star unimpeded. Their only memory of thermodynamics is the thermalized matter that 
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(C.l) 
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produced them. In analogy with the photon entropy production (see Appendix A), the 
neutrino entropy production may be written: 



„ = JdV JdE^ , (C.2) 



where n E is the neutrino differential production rate density and we integrate over the 
single neutrino energy E, instead of frequency. Unlike photons, there is no matter heat 
bath term for free-streaming neutrinos. Technically, there is a second term in (C.2) giving 
the tiny contribution to T, v from the neutrino absorption (as in terrestrial detectors). The 
overall structure of Y, u is non-local because neutrinos are generally not in LTE. 

Neutrinos have an energy- dependent brightness temperature T E , but, if unconfined, no 
chemical potential. This temperature is defined by the Fermi-Dirac analogue of (A. 8): 



■l = ±m 

T E E 



2{E/hcfc 
8tt 3 / r 



- 1 



(C.3) 



where I E = E N E is the energy-specific neutrino differential energy flux; N E is the same for 
neutrino number. Define additional functions: 



h u = J dE h E , N u = j dV n 
N E = JdV h E , N E = N E /4nr 



(C.4) 



where the last definition holds only in the case of spherical symmetry. The exact result 
of (C.4) depends on the specific reaction. (Neutrino entropy production will be further 
discussed in future publications; see Essex and Kennedy 1996.) 



2. Examples 

For an infinitely sharp line, unbroadened by thermal fluctuations, the entropy 
production rate is zero: the neutrinos are produced at exactly one energy. The effect of 
thermal broadening for neutrino line emission can be treated in the same way as it is for 
photon emission lines (Stix 1989). Define the thermally broadened shape function 4> E (E) : 

where H(b, z) is the Voigt function: 

H(b, z) ee (fo/Tr) JZo dy eM-V 2 )/[(z - V? + b 2 } , (c fi , 

JZ dzH(b,z) = ^ ■ 1 j 
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I E and A^ are proportional to <f> E . The reduced variable z = (E — E )/AE D , where E is 
the center of the neutrino line and AE D is the Doppler width: 



AE 



D 



En 



2T 



(C.7) 



V m^C cr 

A is the nucleus emitting the neutrino, and £ t is the root-mean-square microturbulence 
speed, set to zero here for simplicity. Define 7 as twice the effective collision rate for the 
nucleus A in the plasma; then: 

b^T^r- • (C8) 

The collision rate is dominated by collisions of A with the plasma electrons, and, usually, 
6 < 1. If b « 1, H(b, 0) = 1 + (9(6). 



In the spherical case, 

where the background functions are distinguished again by *. 
since we lack neutrino LTE: 

N E = J dV n E {r') . 
A crude estimate of T, u for a broadened line is: 

-(E /hcfR 2 c v D 



(C.9) 



This expression is non-local, 



N v In 



where vp is the Doppler root-mean-square speed: vd = cAE D /E , and R c is the radius of 
the neutrino-producing stellar core. In the Sun, for example, the neutrino lines produced 
by 7 Be contribute approximately 5xl0 18 erg/°K sec to H v . 

If we let T — > holding all other state variables fixed, the entropy production T, u of a 
line vanishes as T^lnT 1 / 2 , as it should: at zero temperature, the reactant matter is not 
thermalized. 

In the continuum case, we distinguish the actual neutrino energy E from the 
zero-temperature neutrino energy E , which is now continuous. We then define a double 
differential neutrino production rate density riEE , per unit dE and per unit dE . Then: 

-(E/hcfr 2 

- 7Y 2 hN EEo 

In the case of the Sun's dominant pp-produced neutrinos, we estimate a contribution to T, u 
of 10 23 erg/°K sec. Because the Sun's photon luminosity so much dominates its neutrino 
luminosity, the solar neutrino entropy production rate is far smaller than that of solar 
photons. 



£„ = JdV J dE J dE [tieeq]* In 



(C.ll) 
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Fig. 1.— Entropy production rate £ of a potentially convective cell conducting heat 
only (solid line) versus conducting heat and rising by convection together (dashed line). 
Convective solution is preferred (lower S) if Schwarzschild criterion is satisfied: V > V a d- 



